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' We present a straightforward, divergence-free derivation of the equation of motion for a 

small but finite object in an arbitrary background using the strong field point particle limit. 
The resulting equation is that of a generalized geodesic for a non-rotating spherical object. 
This is consisitent with the results of previous studies. 
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. Equations of motion in general relativity have attracted much interest recently 

because of the belief that it may be possible to use gravitational waves as an astro- 
' nomical tool in the near future. In fact, several gravitational wave detectors are now 

. operating around the world, ^^"^-^ and there are plans to build the next generation 

Q I detectors. Gravitational waves may open a new window to the universe, just as radio 



waves and X-rays did. Gravitational waves should allow us to study entirely new 
Q^. regimes in the universe, such as the surfaces of black holes and the very beginning of 

^ I the universe, that no other type of radiation can reveal, and they will probably allow 

W)' us to find new astronomical objects that we have not imagined to exist. However, in 

order to realize this, not only technology but also theory must be developed further. 
Theoretically, it is necessary to have accurate predictions of the nature of the waves 



>< 



■ emitted by various astrophysical sources. For this purpose, we need to have a good 



understanding of equations of motion with radiation reaction in a wide variety of 
situations. In particular, at present we have only a very poor understanding of the 
equation of motion for a fast moving source, such as small but extended objects (e.g., 
neutron stars or black holes), in an arbitrary external field. In such a situation, the 
object does not necessarily move along a geodesic of the external field, because the 
self-field of the object and the gravitational waves generated by the orbital motion 
are not negligible. 

In this paper, we study the equation of motion describing such a situation. 
It should be mentioned that Mino et al. and others have derived an equation of 
motion for a point particle with mass m that is represented by a delta function 
source in an arbitrary background, and the equation is interpreted as a geodesic 
equation on the geometry determined by the external field and the so-called tail part 
of the self field of the particle to first order in rn.^^'^-*'^^-* Furthermore, Mino et al. 
used another approach, the matched asymptotic expansion, to obtain an equation of 



typeset using VTpJ^.ds (Ver.0.9) 
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motion without employing the concept of a point particle and thus avoided divergence 
in their derivation. 

Our approach is different from the previously used approaches for the reason 
that wc take the extended nature of the object into account and thus avoid using 
a singular source. This is done by useing the point particle limit developed by one 
of the present anthers. We believe that our approach simplifies the proof given by 
Mino et al. (employing the matched asymptotic expansion) that the Mino-Sasaki- 
Tanaka equation of motion is applicable to a nonsingular source. 

Let us start by explaining the situation we consider and introduce the strong 
field point particle limit, which plays an essential role in our formulation. A small 
compact object with mass m moves with arbitrary speed near a massive body with 
mass M. We assume that the object is stationary, except for higher-order tidal 
effects, so that we can safely ignore emission of gravitational waves from the object 
itself. However, of course, we cannot ignore the gravitational waves emitted by the 
orbital motion of the object. Wc denote the world line of the center of mass of the 
object by z{t) and define the body zone of the object as follows. We imagine a 
spherical region around z{t) whose radius scales as e. Also, we stipulate that the 
linear dimension of the object scale as e^, so that the boundary of the body zone 
is located at the far zone of the object. Thus, we are able to carry out a multipole 
expansion of the field generated by the object at the surface of the body zone. We also 
implicitly assume that the mass of the object scales as e^, so that the compactness 
of the object remains constant in the point particle limit e — ^ 0. This is why we 
call this limit the strong field point particle limit. Then we calculate the metric 
perturbation induced by the small object in this limit. The smallness parameter 
has the dimension of length, and this length characterizes the smallness of the 
object. One may regard it as the ratio of the physical scale of the object and the 
characteristic scale of the background curvature. We assume that the background 
metric g^i, satisfies the Einstein equations in vacuum. Therefore, the Ricci tensor of 
the background vanishes. Since we have assumed that mass scale of the small object 
is much smaller than the scale of gravitational field of the background geometry, we 
approximate the metric perturbation by the linear perturbation of the small particle, 



where the semicolon represents the covariant derivative with respect to the back- 
ground metric, and the trace-reversed variable is defined as usual: 



We employ the harmonic gauge, in which we have 




(1) 



(2) 



Then the linearized Einstein equations take the form 



^h^^'f^ (x) - R"^"^ {x) Up {x) = SttTI"' (x) . 



(3) 
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This can be solved formally as 

h^^'' {x) = Stt y d^y^G^\^ {x, y) T-^ (y) , (4) 

where we have used the retarded tensor Green function defined by 

G'^'^"^ (x, y) = {E {x) , y) \u^'^^P (x, y) 5 (a (x, y)) + i;'^^"^ (x, y) 9 {-a {x, y)) . 

47r L J 

(5) 

For a general tensor Green function and the definitions of a (x, y), y'^" (x, y),u^"''^'^ (x,y), 
y-^^ jj j-^^^ please refer to Refs. ^'and ^' . 
Now we take the point particle limit. In this limit, the above metric perturbation 
contains terms with different e dependences. This simplifies the calculation of the 
equation of motion. For example, terms with negative powers of e appear from the 
delta function part of the Green function, 

,/i^-(x) = 2 J d''y^gu>'\p{x,y)5{a{x,y))T^P [y) . (6) 

As explained below, we only need the field on the boundary of the body zone, 
which is the far zone of the body itself, and thus we may make use of the multipole 
expansion for the field. For this purpose, we choose our coordinates as 

y« = z° (r) + 5y", 

z°(r)=r, <5y° = 0, 

where (r) is the world line of the center of the object. The center is assumed 
always to be inside the body in the point particle limit, and thus there is no ambiguity 
regarding the choice of the center. In these coordinates, the volume element satisfies 
d^y = drd^Sy. Thus we have 

sh^"" (x) = 2 / dHy^gu>'\p (x, z (r,) + 5y) -^T^P (y) , (7) 

where Ty is the retarded time of each point y. Then the multipole expansion is 
obtained by expanding the above expression at the retarded time of the center of the 
object Tz, defined by cr (.x, z (r^)) = 0. This can be done easily by using the condition 
a (x, z (Ty) + 6y) = 0. Then the diference between and Ty is given by 

6r = - '--^'''^f'f +0{5y^). (8) 
a[x,z{Tz)) 

Using this 6t, we can expand u^^ (x, z (vy) + 6y) about r^. In principle, we can 
calculate arbitrary higher multipole moments in this way. Here we only calculate 
the first two terms, namely the monopole and spin terms. Then we only need the 
following expressions: 

"'"'a/3 ix,z{Ty)+Sy) 



a, 
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and 

a {x, z (Ty) + 5y) 
= a{x,z{Tz)) 



J _^ I 0-;7 {x,z{t^)) _ (J(x,z(tJ)0-;^ (x,2:(r^)) ^ 

(t(x,z(t2)) 0-2 (x,2;(t^)) 



(10) 



By defining the mass and the spin as 

mi" (r,) (r,) = j dHy^ T'^^ (y) , (11) 

mS^(" (r,) i^) (t.) = y" dHyV^ Sy^^ T'^^ (y) , (12) 

(see Ref. for the definition of the spin tensor), we finally obtain the following 
expression for gh^'^: 

1m 

O l^x, Z {Tz)) 

2ma{x,z{Tz))a.j{x,z{Tz)) ( WQia/ \-B( \ 

^ ^3 (x,^(r,)) ^'^ a/3(^'^(^-))'5^ (^.)^''(^.) 

Now we derive the equation of motion using this expression. First, we define the 
e dependent 4-momcntiim of the object as the volume integral of the effective stress 
energy tensor Q^^^ over the body zone B{t): 

pt^lj) = - f 0^"'dK- (14) 

JB{r) 

Here we choose the Landau-Lifshitz form for the stress energy tensor, 

0f^^ = {-g)(T^^'^+t1l), (15) 

where T'^'^ is the stress-energy tensor of matter and tj^j^ the Landau-Lifshitz(LL) 
psuedo-tensor, whose explicit expression can be found in their textbook.^'^^ Because 
the effective stress energy tensor satisfies the conservation law O^'^u = 0, the change 
of the 4-momentum defined above can be expressed as a surface integral over the 
boundary of the body zone dB: 



IF 



I e^n^{l + ea-n){-g)t^^j^dQ. (16) 

JdB 



Here, is the unit normal vector to the surface and = Du^/dr the 4-acceleration. 
The equation of motion is obtained by taking the point particle limit, e — >^ 0, in this 
expression. Thus, we need to calculate {—g)t'^L body zone boundary (field 
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points), on which the multipole expansion of the self field and the Taylor expansion 
of the nonsingular part of the field are known, and to retain only terms of order 
in the expression. All the other terms vanish in the limit or the angular integration. 

For simplicity, we demonstrate the derivation of the equation of motion in the 
case of a non-rotating, spherical object. Recalling that the LL tensor is bilinear in 
the Christoffel symbol and that the Christoffel symbol is a derivative of the metric 
tensor, it is found that the only remaining terms come from the combination of the 
0-th order of the smooth part of the metric and the part of the self field proportional 
to e~^, which takes the following form: 

^-Th^'ix) = -7:;^^^^^;3(^,^(r))i'^(r)i^(r)|,=,,. (17) 

(J[X, Z[T )) 

The remaining part of the self field is the so-called tail part, which is regular in the 
limiting process. The field point (x, Tx) is now on the body zone boundary, which is 
defined by ^/2a{x, z{tx)) = e and 

[a.a{x,z{T))z^{T)]r=rx=0. (18) 

We have assumed a stationary spherical source, and in this case, m is the Arnowitt- 
Deser-Misncr (ADM) mass that the compact object would have if it were isolated, 

m = limP'^. (19) 

Using the above expressions in the LL tensor, we have the following result in the 
point particle limit: 



mrUgs)u''u(' - m-FMu^u^^- (20) 



Here the Christoffel symbol is evaluated using the 0-th order of the smooth part of 
the metric Qs- Using the fact that the ADM mass is related to the 4-momentum as 



= V^smui" (21) 

(which is suggested by the higher-order post Newton approximation-*^^)), we finally 
have 

^ = -r'^,{9s)u-u^ (22) 

which is the geodesic equation on the geometry determined by the smooth part of 
the metric around the compact object. 

In fact, the spin effect on the equation of motion can be derived in a similar way, 
and the standard result is obtained. 

We have proved that a small compact object moves on the geodesic determined 
by the smooth part of the geometry around an object. The equation describing this 
motion is the so-called MiSaTaQuWa equation. The smooth part contains gravita- 
tional waves emitted by the orbital motion, and therefore the equation contains the 
damping force due to the radiation reaction. It is necessary to explicitly calculate the 
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smooth part of the metric generated by the small object around the orbit. This is the 
heart of the problem, and it is not attempted in this paper. However, we point out 
that our method avoids use of a singular source, employing only the retarded Green 
function, by making use of the point particle limit, and all the quantities are evalu- 
ated at the surface of the body zone boundary. Thus we only need the dependence 
of the distance from the center of the object, namely, the e dependence of the field. 
For these reasons we arc able to avoid using any divergent quantities in any part of 
our calculation. One of the difficulties encountered when using a singular source is 
the ambiguity of the separation between the divergent and convergent parts at the 
singularity. With our method, it might be possible to avoid this kind of ambiguity 
and to obtain a unique equation for the fast motion with a radiation reaction. We 
would like to see this if this is indeed the case in a future work. 
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